Theoretical study of the quenching of NH (^A) molecules via collisions with Rb 

atoms: preliminary results 
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We examine the quenching reaction Rb (^S) + NH (^A) Rb (^Pi/2) + NH {X ^E"). This 
reaction may be utilized to produce ground state NH molecules for studies of ultracold physics or for 
other purposes, and is interesting in that it involves initial and final states that are nearly degenerate. 
This near degeneracy is expected to lead to a large reaction rate. We examine this system using ab 
initio quantum chemistry calculations and scattering calculations, which include spin-orbit effects, 
and find that the reaction rate is large and, in fact, approaches the quantum mechanical unitarity 
limit. We discuss the prospects for an experimental examination of this system. 

PACS numbers: 



I. INTRODUCTION 

Chemical reactions involving free-radical molecules are 
important in many dense gas systems including combus- 
tion gases [ij and interstellar clouds 0, 01 . The fundamen- 
tal imidogen free radical(NH) has been detected [J as an 
intermediate in the combustion of CH4 with N2O. NH has 
a ground state triplet, X^Yj~ , and a long lived metastable 
stateQ, the singlet delta, labeled aA^ or ^A. This state 
is doubly forbidden to decay via electric dipole radiation 
(^A X^Y,^). Therefore, it is a good system in which 
to study properties of excited molecular states @. 

Currently, there is considerable interest in collisions be- 
tween diatomic molecules, such as NH, and alkali atoms. 
This interest stems from the relevance of such systems to 
ultracold physics and chemistry 0, H, H, [13, [ill, [H, [11 
and the increased level of control possible through the use 
of cooling and trapping techniques. The advent [l^ of 
the Bose-Einstcn condensate (BEC) led to the produc- 
tion of mixed BECs [l^ and engendered work on pro- 
duction and trapping of cold molecules [H, [13, [H, O ■ 
There has been considerable recent experimental work 
on creating cold dipolar molecular samples for colli- 
sion studies. For example, cold OH and NH 
molecules in either the ground ('^S~) or the ^A elec- 
tronic state [21I, [23, m, [Ij, [2^ have been produced. In 
addition, there are a growing number of theoretical in- 
vestigations on interactions between diatomic molecules 
and alkali atoms. A broad survey of interactions of al- 
kali and alkaline earth metals with NH was published 
in Ref. [2^. Other theoretical work on the Rb-NH sys- 
tem involved examinations of the relevant potential en- 
ergy surfaces [13, [H, [11]. Explorations have also been 
conducted for related systems such as Rb-OH Isol Isil . 
Hc-NH [3i,[3l,[33 Rb-NHg and NH-NlTjU"^ 

Our interest aims toward a broad treatment of the Rb- 



NH system, including up to the first excited electronic 
states of both the atom and molecule. We have con- 
structed a preliminary treatment and present some ini- 
tial results here. The goal is to accurately calculate the 
dynamics on the higher excited electronic states. In par- 
ticular, we concentrate on the states ^^S^ and ^A of NH 
and the states ^S", ^Pi/2, and ^^3/2 of Rb. As it hap- 
pens, the excitation energies from the state to the 
states are nearly equal to the excitation energy of NH, 
such that the channels are unusually close in energy; in 
fact, the excited NH channel lies between the two spin- 
orbit components of the Rb ■^P state. Specifically, the 
excitation energies are 12687. 8cm^^ on the NH [3^, and 
12578.96 and 12816. SScm^i for Rb ^P^/a and ^Pg/j, re- 
spectively [soj . 

This situation makes the quenching of the excited NH 
(^A) state by Rb, using the reaction 



Rb(2S) + NH(iA) 



Rb(2p) + NH(X^S") , (1) 
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an interesting one to study for several reasons. First, 
it opens up the possibility that this quenching reaction 
could find some use in ultracold physics, as it would rep- 
resent an electronically-inelastic reaction at cold initial 
and final collision energies. This is a chemically-reactive 
process, in the sense that the chemical nature of the prod- 
ucts is dramatically different from that of the reactants. 
Second, it could prove useful as a mechanism for produc- 
ing ground state NH from the exited state. 

This quenching reaction can be precisely studied us- 
ing techniques to cool and trap atomic and molecular 
samples. The methods of atom cooling and trapping 
have been refined over the last 15 years such that one 
can routinely produce dense samples (up to 10^'^ cm~'^) 
with temperatures on the order of 1-100 ^K. More re- 
cently, with the development of Stark deceleration of po- 
lar molecules [40j. trapped molecular samples can be pro- 
duced with temperatures on the order of 1 -100 mK. 

An experimental realization will proceed as follows. 
The atoms will be cooled in a magneto-optical trap and 
transferred to quadrupole magnetic trap in one region of 
the vacuum system. During this process, a pulsed beam 
of NH molecules will be slowed using a Stark decelera- 
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tor and trapped using electrostatic fields [4l| . Because of 
the weak interaction between the atom and the electric 
field, the two species can be controlled independently. 
This control will allow the atoms to be moved to over- 
lap the molecular sample, setting the initial time for the 
interaction. There will be several methods of detection 
employed in order to gain a full understanding of the 
interaction. A novel signature of this reaction is the pho- 
ton produced from the decay of the excited atom (^Pi/2 
—y ^81/2)- This dipole-allowed transition will happen 
very rapidly after the coUisional excitation of Rb. These 
photons can be detected efficiently with almost no back- 
ground using a spectrally filtered photomultiplier tube. 
In addition, both the number and temperature of the 
trapped atom sample can be measured with absorption 
imaging; the number and temperature of the molecular 
sample, with REMPI. Through this procedure, it should 
be possible to measure the reaction rate. 



II. QUANTUM CHEMISTRY CALCULATIONS 

We carry out quantum chemistry calculations using 
the COLUMBUS package [H, [H, 13] for calculations on 
the electronic states of the Rb-NH system. We use an 
effective core potential [ill, that accounts for the 36 
core electrons of Rb, and thus treats the Rb atom as a 
one-electron system. We therefore perform calculations 
on nine-electron states. We adjust our effective core po- 
tential to reproduce the correct channel energies of the 
combined Rb-NH system at infinite separation. 

For the moment we have constructed a set of two- 
dimensional potential energy surfaces, fixing the NH 
bond length at 1.925ao- We refer to the configuration 
of the Rb-NH molecule using Jacobi coordinates. These 
coordinates are defined as r, the NH bond distance, 
1.925ao; R, the distance between the NH center of mass 
and the Rb; and 7, the angle between the two corre- 
sponding vectors, such that 7=0 denotes a linear Rb-H-N 
configuration. 

Owing to our use of a pseudopotential, in addition to 
the large number of electronic states treated by a rela- 
tively modest configuration-interaction (CI) calculation, 
the surfaces we construct are not expected to be of "spec- 
troscopic accuracy" but instead probably have errors on 
the order of tens of meV. They are expected to be suffi- 
ciently accurate for the qualitative study presented here, 
and in any case, the omission of motion in the Jacobi rjyH 
degree of freedom surely compromises the calculated dy- 
namics more than does any error in the potential energy 
surfaces. 

In all, and without considering spin-orbit coupling, 
there are eight doublet states and four quartet states, 
for a total of 32 microstates. This odd electron system 
enjoys the Kramers degeneracy and thus the electronic 
Hilbert space is block-diagonal in two 16x16 blocks. For 
total angular momentum J ~ we need only consider 
one of these blocks. The asymptotes of the adiabatic 



TABLE L Electronic microstates considered in our study: 
quantum numbers at linear geometry. Each is a member of a 
different Kramers doublet. The quantum numbers also label 
the members of our diabatic basis for all geometries and the 
electronic channels of the scattering calculation. The quan- 
tum numbers \Q\ and S are the absolute values of the projec- 
tion of the total and spin angular momentum on the NH axis. 
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states are listed in Table HI 

At linear geometry, the electronic states are described 
by their projections of angular momentum on the molec- 
ular axis. We refer to the electronic states by the stan- 
dard A, 17, and S quantum numbers. A is the absolute 
value of the projection of electronic angular momentum 
on the molecular axis, the eigenvalue of the operator Iz, 
such that A=2 denotes a A state, etc. \n\ is the absolute 
value of the projection of the total angular momentum 
on the molecular axis, the eigenvalue of the operator jz 
= Iz + si. The absolute value of the eigenvalue of si 
is denoted E. The symbol E is also used to denote a 
state with A=0, but we will avoid ambiguity by always 
denoting the multiplicity of a given electronic state, as in 

2E. 

The Born-Oppenheimer electronic states are first cal- 
culated without spin orbit coupling. The next step is a 
diabatization [47l | of these potentials, in which we con- 
struct diabatic states labeled primarily by integer values 
of A, by diagonalizing the electronic angular momentum 
Iz projected onto the NH axis. ±A is thus not only the 
eigenvalue of the projection of angular momentum at lin- 
ear geometry for a given diabatic electronic state, but 
also its approximate expectation value at other geome- 
tries. We have a pair of diabatic states, two pairs of 
^H states, and two states. The doublet H states are 
distinguished as anionic in character (NH~ x Rb+ ) or 
not (NH X Rb 2p). The E states are essentially NH 
■^E^ times Rb in either the ground state or the ^P^ 
state. We also have two *E and one ''H state. 

A diabatization is useful to us for two reasons. First, it 



3 



allows for the convenient inclusion of nonadiabatic effects 
among the adiabatic surfaces, accounting for conical in- 
tersections, for instance. Rovibronic effects such as the 
Renner-Teller effect can also be accounted for by such a 
diabatization, though at this stage we do not include the 
mixing of electronic and rotational angular momentum. 

Second, it appears to allow us to add the spin-orbit 
terms "by hand," and therefore enables us to perform a 
larger electronic structure calculation, though we employ 
a smaller calculation using the spin-orbit CI capability 
of COLUMBUS to verify the accuracy of this procedure. 
The final steps in constructing the surfaces are then the 
addition of spin-orbit terms to the diabatic electronic 
Hamiltonian and, last, transformation to a representa- 
tion that accounts for the Kramers degeneracy. 



A. Details 

We have found it prohibitive to use a basis set large 
enough to get a sufficiently accurate NH excitation en- 
ergy. Therefore, in order to reproduce the physical ener- 
getics, we artificially modify the effective Rb core poten- 
tial. We beg in with the SDF pseudopotential developed 
in Refs. [45l. l46j. The pseudopotential for the S and P 
waves is then modified by changing their functional forms 
into: 



45.272 exp(-1.012r2) ^ 45.272 exp(-1.14973057r2) 
P ■■ 

2.83cxp(-0.3036r2) 2.83 exp(-0.267230608r2) 

(2) 

An uncontracted SDF basis set is adopted for the Rb; 
while for N and H we use the aug-cc-pvtz basis set of 
Dunning [i^. The excitation energies obtained are ap- 
proximately 1298(5) cm~^ on both the Rb (without spin- 
orbit interaction) and NH fragments. 

The electronic configuration of the NH states of inter- 
est is ls^2s^2p^2pj. In addition to these orbitals we have 
the four orbitals on rubidium (s and p) , for a total of nine 
orbitals in the valence space. 

To describe all of these states for all geometries, it is 
necessary to include orbitals beyond the minimum set of 
nine, 1 — 5(t and 1 — 27r, orbitals. The anion state is the 
highest in energy at infinite Rb-NH separation, but be- 
comes the ground state at small Rb-NH separation, and 
it undergoes avoided crossings with all of the other states, 
and so is relevant to the dynamics on the other surfaces. 
To describe the anion state, two additional orbitals are 
necessary: a tt relaxation orbital that accounts for the 
expansion of the NH tt orbitals in the anion state rela- 
tive to the neutral states; and a correlating a* orbital on 
the NH which increases the Rb-NH bonding on the an- 
ion surface. In Ref. [13] a similar orbital space was used 
although what we denote as an anion relaxation orbital 
was denoted by those authors as a 6p orbital. 



For the first step, a state-averaged multiconfigura- 
tion SCF (MCSCF) calculation in the minimum 9-orbital 
space is performed on the lowest six doublet states and 
the four quartets. This calculation produces a tt orbital 
on the NH fragment that docs not account for the relax- 
ation within the anion state; that calculation significantly 
overestimates the anion state energy at small Rb-NH sep- 
arations. 

The next step is an all-singles-and-doubles 
configuration-interaction (CI) calculation on the eight 
doublets and four quartets in which we are interested. 
At large R, we must follow the anion state surface as it 
rises above several states that we are not interested in. 
Twelve averaged natural orbitals are adopted from this 
CI calculation, expanding the orbital space to include 
both the TT relaxation orbitals and the a* correlating 
orbital. We then perform an all-singles-and-doubles 
CI calculation in this space of twelve natural orbitals. 
In both CI steps we freeze the l-2a', a.k.a. the l-2cr, 
orbitals on NH, which are basically the N Is and 2s. 
Among the four sets of spin and space symmetries the 
maximum number of configuration state functions in the 
CI is 1,110,312 for the doublet A'. 



B. Diabatization 

At this stage of the calculation there are twelve adi- 
abatic electronic states as functions of nuclear geome- 
try (eight doublets and four quartets), and the spin-orbit 
terms have not yet been added. A unitary transformation 
is now applied to these twelve states, a property-based 
diabatization [13, [H, [l^l ■ 

The doublets are diabatized using both the electronic 
dipole operator and the z-projection of the orbital an- 
gular momentum relative to the NH bond axis (the 1^ 
operator with eigenvalue ± A at linear geometry) . The 1^ 
operator is first diagonalized in the adiabatic basis. In 
order to prevent an avoided crossing at large Rb-NH sep- 
arations that would lead to a problematic mixing of our 
diabatic basis, a column and row of the Iz matrix must be 
damped. This transformation of the Iz operator does not 
affect the eigenvectors, except near the avoided crossings. 
The matrix elements of Iz in the ground-state A" row and 
column are thus reduced by hand beyond R~lOAaQ. At 
such large internuclear separations the lowest A" state is 
Rb (^S) + NH (^S~). We reduce the matrix elements in 
this row and column by a geometry-dependent factor 

(C),i--(C).i fl i?sin^^ < 10.4ao 

V, 1 1 X <^ 10.4 (3) 

{lz)i,t ^ {lz)i,t yRsinO i?sin0 > 10.4ao . 

With this damping function, the mixing between the 
two states NH (^S") times Rb {^S) or (^Pa") is sup- 
pressed and the nonionic H state has a Iz eigenvalue 
which plateaus around 1.85. 
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FIG. 1: Eigenvalues and expectation values of the operator h 
for diabatic states evaluated at i?_Ri,_jv_r/ = 16.0ao as a function 
of the Jacobi angle 7. (a) Eigenvalues in the adiabatic basis, 
(b) Eigenvalues in the adiabatic basis with damping function 
of Eq|31 (c) Expectation values of Iz for final diabatic states. 
The absolute values of these numbers provide the A labels of 
the diabatic states (^E, ^11, ^A). 
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FIG. 2: a) Non-spin-orbit CI results at 7 = 180°. b) Final 
spin-orbit surfaces used in study, obtained by adding the spin- 
orbit term by hand to the diabatic basis, 7=180°. 



similar to that of RbOH [30|, l3l| . Diagonalization of the 
electronic dipolc operator in the direction of the Jacobi 
vector i? is carried out in the 11+ and spaces sepa- 
rately, and the electronic Hamiltonian is diagonalized in 
the space. The quartets are diabatized by diagonal- 
ization of Iz , followed by diagonalization of the electronic 
Hamiltonian in the "^S space. The expectation values of 
Iz in our final diabatic basis are all near integers. 



Fig.[T]shows the eigenvalues of 1^ for the doublets only. 
FigHh depicts the eigenvalues of Iz in the adiabatic ba- 
sis at i?=16.0ao. The Iz eigenvalues using the damping 
function are plotted in Fig[T]3. The final expectation val- 
ues of Iz in our diabatic basis are shown in Fig. [T]:; these 
values are nearly integers, which permits the labeling of 
these states by the A quantum number labels. We have 
a total of two ^S, two ^11, and one ^A, and two "'E and 
one *n diabatic state. 

The electronic dipole operator is used to separate the 
two states. These correlate with the anion NH~ (^11) 
+ Rb+ state and with NH [X^Y,-) -t- Rb (^P) states at 
infinity; at small Rb-NH separations the former is the 
lower electronic state and at large separations the lat- 
ter is. The interaction between these states is such that 
they form an avoided crossing of width approximately 
1000 cm""'^ around Rm-NH — llao- The situation is 



C. Adding spin-orbit terms 

At this point, then, we have a diabatic basis in which 
each member is labeled by a given projection of elec- 
tronic angular momentum about the NH bond axis, and 
in which the spin-orbit part of the electronic Hamilto- 
nian has not been included. Multiplying each of these 
electronic basis functions by all possible spin wavefunc- 
tions, we arrive at the full 32-microstate basis. Assuming 
that our diabatic basis function labels are indeed good 
quantum numbers, the addition of the spin-orbit terms 
is straightforward, requiring only the enforcement of the 
Condon-Shortley phase convention among the members 
of the diabatic basis, and the spin-orbit terms have no 
dependence upon the nuclear geometry. The discussion 
below verifies that this approximation faithfully repro- 
duces the true eigenvalues. 

After adding these spin-orbit terms we shift the delta 
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Jacobi Rpb-NH (units of bohr) 

FIG. 3: Long-range behavior of surfaces, 7 = 180°. 

diabatic potential energy surfaces slightly (by tens of 
wavenumbers) such that the asymptotes coincide with 
their physical values. A unitary transformation now 
yields the two uncoupled 16x16 blocks each containing 
one member of every Kramers doublet. These blocks are 
complex conjugates of one another; it is only necessary to 
choose one of these to use as the Born-Oppenheimer elec- 
tronic Hamiltonian for total angular momentum J = 0. 
Extensions to include the electronic angular momentum 
(the Renner- Teller and spin Renner- Teller effect) and 
nonzero J will instead use the real- valued 32x32 rep- 
resentation. 



D. Description of the surfaces 

Results of the non-spin-orbit CI for linear Rb-N-H ge- 
ometry are shown in Fig. [2^. The ground state in the 
asymptotic region is Rb (^S) + NH (X^S^) with an 
asymptote of -55.1993162 hartree. 

The anion surface comes sweeping down, with a 
minimum at 5.1192ao of -55.2346523 hartree, giving a 
depth of 0.949eV. In comparison, Ref. [131 found a well 
depth of 1.372eV at i?=4.911ao, for a NH bond length 
of 1.948ao. It crosses the collection of nearly degenerate 
channels at approximately 12.0ao, and crosses the sur- 
face that is asymptotically the ground state at 7.3029ao 
and at an energy of -55.1997712 hartree. In comparison, 
Ref. [231 found the crossing at 7.163ao, also for a NH 
bond length of 1.948ao. 

The final results including the spin-orbit interaction, 
obtained by adding the spin-orbit terms by hand to the 
diabatic Hamiltonian, are plotted in Fig. [2j3, in a closer 
view. The state is visible as the flattest line, staying 
near zero wavenumbers until around lOao on the scale 
of the figure. At the far edge of the figure the surfaces 
have separated into the Rb (^S) x NH (^A) asymptote, 
in the middle, with Rb (^P3/2) and Rb CP1/2) times NH 
(X^E") above and below, respectively. 

Components of the upper ^P3/2 state split into sets 




10 15 20 25 30 
^Rb-NH (units of bohr) 



FIG. 4: Spin-orbit results using a smaller, 9 orbital va- 
lence space CI, 7=30°. Spline-interpolated surfaces calcu- 
lated from diabatic states adding the spin-orbit interaction 
by hand (lines) compared with spin-orbit CI results (points). 
The asymptotes are not correct for this test calculation. 



following either the ^'^E surface, which rises from its 
asymptote inward, or the ^'^H surfaces. These latter sur- 
faces drop in energy going inward from their asymptote 
- the quartet less so, whereas the ^H surface undergoes 
an avoided crossing with the anion ^H surface, dropping 
steeply downward from the top of the figure. Inward of 
this avoided crossing around 12.0ao, the ^H (Rb x 
NH X^'S^) state is found dropping downward with de- 
creasing R from the avoided crossing and crosses the ^A 
surface around 9.0ao- There are other crossings farther 
in, too numerous to mention. 

The long range behavior of the surfaces is shown in 
Fig. [31 This behavior is governed by the long-range in- 
teractions: as the NH has a dipole moment, these are 
dipole-quadrupole with ^ power law for the Rb (^P) 
fragment and dipole-induced dipole with ^ power law 
for the fragment. The dipole-quadrupole interaction, 
acting on the Rb ^P x NH '^E~ channels, is stronger than 
the dipole-induced dipole interaction for the ^A channel 
and therefore the crossings are determined by the for- 
mer. At linear Rb-N-H geometry, this interaction splits 
the ^P state into the sigma component at higher energy 
and the pi component at lower energy, which behavior we 
described above for the results in Fig. ^p. It is responsi- 
ble for the crossings around 14.7ao. We examined these 
crossings for various Jacobi angles 7; for all angles the 
crossing is apparently too sharp to provide a mechanism 
for coupling between the ^A and the components of the 
Rb (^P3/2) state that correlate with H symmetry farther 
in. Results below seem to indicate that this is indeed the 
case. The relevant crossings for collisions at low energy 
are those visible in Fig. [^b. 
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E. Verification of treatment of spin-orbit effect 

The spin-orbit calculation using the COLUMBUS pro- 
gram was prohibitively large to perform in the full orbital 
space. We thus perform the spin-orbit calculation only at 
the first, nine-orbital CI step, as described above. Results 
of this calculation are presented in Fig. 21 These results 
verify that addition of the spin-orbit terms by hand to the 
diabatized Hamiltonian from the non-spin-orbit CI cal- 
culation reproduces the results of the full spin-orbit CI. 
(The asymptotes are not correct for the test calculation 
shown in this figure.) 



III. SCATTERING CALCULATION 

We calculate the quantum nuclear dynamics on the 
coupled set of diabatic potential energy surfaces for to- 
tal angular momentum J = 0. The standard [5l|, [s^ 
BF Hamiltonian for R times the wavefunction, keeping r 
fixed, in Jacobi coordinates is 



H 



1 



1 d 



? + U(i?,r,7) 



. sin(7) — 

sm(7) 07 07 



(4) 



where hr is the reduced mass in that degree of free- 
dom; Br is the rotational constant of NH, taken to 
be 16.699cm^^ the value for the ground electronic 
state; and V is the matrix representation of the electronic 
Hamiltonian in the diabatic basis. 

The scattering calculations employ the R-matrix prop- 
agator technique of Baluja, Burke and Morgan [s^l- Our 
implementation adopts the discrete variable representa- 
tion (DVR) [H, [Hi], with the Legendre DVR [Ilin 7 
- with 80 points - and the Gauss-Lobatto DVR [58[ in 
R, with 6 points per element, 960 elements, from 3.25ao 
to 43.25ao. The propagation covers one element at a 
time. For each diabatic electronic state, the basis in 7 
is contracted by calculation of an adiabatic basis in 7 as 
a function of the scattering coordinate R and include 42 
adiabatic states in 7 per diabatic electronic state. Slow 
variable discretization (SVD) [s^ efficiently accounts for 
the coupling between these surfaces nonadiabatic in R. 

Results of the scattering calculation are shown in 
Fig. m The total probabilities for transitions from the 
ground rovibrational state of NH (^A) plus Rb (^S) to 
the other electronic states arc plotted as functions of col- 
lision energy in the incident channel. Two degenerate 



channels in the calculation have Ifil = 



or 



but the 



2 " 2' 

method permits the evaluation of four independent cross 
sections, for both the calculated process and the time- 
reversed process corresponding to the the other 16x16 
Kramers doublet block (i.e., complex conjugating the 
electronic hamiltonian in the diabatic basis). The four 
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-A \a\ = 3/2 TR' 




elastic 



50 100 I 150 200 

Collision energy, wavenumbers 



FIG. 5: Results of scattering calculations: the probability per 
collision for a transition from the lowest rovibrational state of 
(NH X Rb ^S) to other electronic states as a function 
of collision energy, for J = 0. The four unique initial states 
are described in the text. The thick red (dark) line is the 
probability for transition to any of the Rb (^P) channels; the 
thin green (light) line is that for a transition to any of the 
ground ^'^E channels; the dots mark the probability of staying 
in the ^ A channel space. The arrow marks the energy at which 
the ^P3/2 channels bcome open. 



corresponding sets of scattering amplitudes correspond 
to two rows and two columns of the calculated S-matrix. 
The S-matrix cannot be chosen to be symmetric because 
we have excluded the time-reversed orthogonal comple- 
ment of our sixteen states, the sixteen states of the other 
Kramers doublet block. The label "TR" in Fig.Odcnotes 
the time-reversed partner. 

These results indicate that the cross section for the 
quenching reaction is indeed large in the present treat- 
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ment, and it is in fact comparable to the unitarity limit. 
In contrast, the cross section to the ground ^'''E chan- 
nels is significant for low collision energy but is at least 
an order of magnitude below the quenching reaction cross 
section at energies above the Rb ^P3/2 threshold, which 
is marked with an arrow on the abscissa of Fig. [S] A 
qualitative change in the branching ratios becomes evi- 
dent near this energy, 128 wavenumbers; the lower three 
panels of Fig. [5] indicate that for collisions in those inci- 
dent channels, the proportion of Rb (^P) -|- NH (^E~) to 
Rb (^S) + NH (^A) produced increases at higher energy. 
However, the correlation of features such as the change 
in this proportion with the opening of the ^^3/2 channels 
is not perfect, and there is the possibility that it derives 
from the opening of different channels (for instance, the 
Rb (^S) + NH (^A j = 2) channel at 100 wavenumbers, 
or Rb (2Pi/2) + NH (^E" j ^ 3) at 91 wavenumbers), 
or another mechanism altogether. A fuller characteriza- 
tion of these results will be presented in a forthcoming 
publication. 



IV. TIME-DEPENDENT CALCULATIONS 

Time-dependent calculations help to illustrate the dy- 
namics that drive the quenching reaction. These have 
been carried out using a smaller basis of 20 adiabatic 
basis functions in the Jacobi angle 7 using 50 Gauss- 
Legendre DVR basis functions, due to memory con- 
straints, compared with 42 and 80 for the converged R- 
matrix calculation. The propagated incident gaussian 
wavepacket, with a width of \/2 bohr in the R degree of 
freedom, starts in the lowest adiabatic (in R) channel of 
the ^A = I electronic state, at a radius of R=22ao. 
The wavepacket is given an initial translational energy of 
90 cm~^, such that the lower Rb ^Pi/2 channels are open, 
but the ^P3/2 channels are closed, though the bandwidth 
of the initial wavepacket is 18 cm~^. We employ a Lanc- 
zos propagator of order 12. 

The results are shown in Figs. [S] and [71 Fig.[n]shows the 
probability density as a function of time and the scatter- 
ing coordinate R, summed over the adiabatic states (in 
R) for each diabatic electronic channel. In this figure, the 
initial state corresponds to the upper left hand panel. In 
Fig. [71 the initial state is colored green (online). Both 
plots show that the density on the initial state electronic 
surface does not penetrate much beyond i?=10.0ao; in- 
stead, it couples to several of the Rb surfaces and also 
reflects. The cut of the surfaces at linear geometry in 
Fig. [2Jd shows the innermost extent of the attractive part 
of the initial state ^A surface, near R=7.0ao, whereas for 
the opposite Rb-H-N geometry the potential well extends 
only to approximately lO.Oap. It is therefore likely that 
the reflected part of the initial wavepacket, most clearly 
visible in Fig. [SI comes from configurations near linear 
Rb-H-N geometry, and that the remaining flux is lost to 
the other surfaces within lO.Oao- 

Some immediate coupling evidently occurs among the 



initial state and all of the Rb (^P) states, as is clear from 
the similar shape of the lowest contour line in the corrse- 
ponding panels of Fig. [51 around 1000 fs. The couphng 
seems to be strongest between the initial ^A state and 
the Rb (^Pi/2) X NH (^E" E^-^^ = 0), \n\ ^ \ state, the 
eighth state in Tableland Fig. [SI The transition between 
the initial Rb (^S) x NH (^A) and this one is the only 
transition corresponding to a conservation of the projec- 
tion of spin angular momentum on the NH, E^^, along 
with a transition to the Rb ^Pi/2 state. This electronic 
transition does not conserve the projection of the elec- 
tronic angular momentum on the NH axis and therefore 
must be driven at nonlinear geometries. 

Fig.[71shows the density on each adiabatic (in R) curve, 
colored (online) according to the electronic channel in- 
dex. In viewing this figure one should keep in mind that 
there is electronic coupling between these curves, which 
represent the energies of adiabatic (in P) basis functions 
calculated on each diabatic electronic surface. In partic- 
ular, the coupling among some of the excited state curves 
and the ground state sigma curves means that there is 
repulsion on the excited state curves and attraction on 
the ground state curves that is not represented in this 
figure. 

This figure shows that the coupling is indeed strongest, 
earliest, to one of the Rb (^Pi/2) channels, that which 
rises, from large P to small P, around 12ao, and that 
correlates to Rb (^Pi/a) + NH (^E" j 1). Although 
some red and orange (online) are visible on the left hand 
side of this figure, showing the large view including all of 
the corresponding anion and ground state sigma curves, 
the bulk of the density resides clearly on electronic states 
involved in the quenching reaction. 

Despite the fact that the wavefunction amplitude on 
the anion ^H curves never achieves a large value, these 
states have significant influence upon the dynamics, as 
would be expected from the large avoided crossings that 
they create with some of the Rb ^P surfaces in the 
electronically adiabatic picture. R-matrix calculations 
performed without these diabatic states produced elec- 
tronically inelastic cross sections markedly different (and 
smaller) than the ones shown in Fig. [5l The coupling to 
the anion surfaces clearly involves a large amount of rota- 
tional excitation of the NH fragment, as can be seen from 
the fact that all of the adiabatic curves corresponding to 
this electronic state have amplitude on them, colored red 
(online) in Fig. [71 

Because the time-dependent calculation shows that 
the transition is strongest from the initial state to a 
Rb ^Pi/2 state, it corroborates our estimate that the 
avoided crossings around 14.7ao between components of 
the higher-energy, closed Rb ^P3/2 surfaces and the ^A 
states, caused by the long-range dipole-quadrupole in- 
teraction, are sufficiently sharp such that they do not 
provide a good coupling mechanism for the quenching 
reaction (at least for the energy considered in this time- 
dependent calculation). Instead, the evidence suggests 
that the avoided crossings and conical intersections in- 
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FIG. 6: Results of time-dependent wavepacket propagation described in the text. The density, in arbitrary units, integrated 
over 7, on each electronic channel is shown as a function of the time of propagation and the value of R. The results for the 
last state are the same as the one shown, i.e., insufficient density to be visible on the plot. The first contour is at ^^jth the 
value of the second one, and the spacing is linear. The states are in the same order (left to right, top to bottom) as in Tabled 
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FIG. 7: Time-dependent wavepacket propagation described in the text. Each frame is a snapshot of the wavepackets on each 
adiabatic curve. The curves are drawn with black hues and the wavepackets on each of these curves is drawn above and below 
the corresponding curve. Each adiabatic curve corresponds to a different diabatic electronic channel and these are distinguished 
by the color of the wavepacket (online only). Green denotes the channels including the incident one. Blue denotes Rb (^P) 
channels; red denotes the anion state; and orange denotes the ground state and ''E channels. The wavepacket is incident in 
the ground rotational '^A \Q\ =3/2 channel. 
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volving the surface and the Rb surfaces further 
in, caused by the fact that the anion surface drops down 
steeply and undergoes broad avoided crossings with these 
surfaces, drive the wavepacket onto the the Rb sur- 
faces. A more thorough analysis is deferred to a forth- 
coming publication. 

The time dependent treatment corroborates the result 
of the R-matrix calculations, that the coupling to the 
ground NH x Rb (^S) electronic state is relatively 

small, whereby most of the outgoing flux avoids these 
channels. The quartet components of this state seem un- 
involved in the quenching reaction for the energy studied, 
as is clear from Fig. [6] A modest flux is seen in the dou- 
blet ground state sigma component, much less than in the 
Rb channels, in agreement with the results of the R- 
matrix calculations. Thus, it appears that the quenching 
reaction is both efficient and selective. 



V. CONCLUSION 

These preliminary results show that the quenching re- 
action Rb (2S) -h NH ( 1 A ) ^ Rb CP1/2) + NH (X^S") 
does indeed proceed with high probability, as may have 
been expected due to the near degeneracy of the elec- 



tronic states involved. The reaction is selective, yielding 
markedly less ground electronic state products. The cou- 
pling mechanism is seen to involve the complicated set 
of avoided crossings and conical intersections that occur 
inwards of 14.7ao, not those due to the long-range inter- 
actions. A more complete study will be published in the 
future, thoroughly analyzing the dynamics that occurs 
between the many electronic states involved in the sys- 
tem, analyzing the effect of the mixing of electronic and 
spin angular momentum with rotational angular momen- 
tum (the Renner- Teller and spin Renner- Teller effects), 
and providing reaction rates for comparison with exper- 
iment. 
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